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ABSTRACT 

-Sh  e  repeated  passage  of  pressure  waves  past  a  given  point,  following  an  explosion  in  an 
enclosure,  results  in  the  long  and  complicated  pressure  time  histories  called  "complex  waves."  The 
waves  are  not  only  complex  because  of  their  temporal  variation,  but  their  effect  on  biological 
structures  depends  on  the  direction  of  each  wave  component.  The  ability  to  estimate  injury  under 
these  circumstances  represents  one  of  the  greatest  challenges  of  biomechanical  modeling. 

mathematical  model  for  the  reflections  of  blast  waves  within  an  enclosure  has  been 
constructed  based  on  the  concept  of  mirror-like  reflections  of  the  wave  front  at  each  solid  surface. 
The  model  accounts  for  the  position  and  orientation  of  each  surface  and  allows  for  the  presence  of 
openings,  such  as  doors  and  windows.  A  simplified  model  for  the  venting  of  explosion  product 
gases  is  also  included*  f  <■  vl  vKicrrVjy  J  ►J  it  ^ 'j  - 

The  pressure  loading  time  histories  are  combined  with  simplified  models  of  the  structural 
response  of  various  organ  systems  that  are  known  to  be  injured  by  air  blast.  The  combined  models 
not  only  allow  estimates  of  potential  injury  due  to  explosions  within  enclosures,  but  also  provide  a 
means  of  judging  the  accuracy  of  the  blast  calculations  in  terms  of  the  requirements  of  injury 
prediction. 

The  blast  field  predictions  are  compared  with  measurements  made  inside  an  Armored 
Personnel  Carrier  and  inside  a  field  bunker.  Many  qualitative  features  are  captured  but  the 
agreement  does  not  extent  to  individual  details.  A  comparison  of  the  peak  overpressures  and  the 
total  pressure  impulse  shows  that  in  most  cases  agreement  is  within  20-30%,  but  there  are  certain 
cases  for  which  the  disagreement  is  much  greater. 

The  measured  and  calculated  pressure  histories  are  used  in  the  simplified  biological 
response  models  and  a  comparison  made  of  the  critical  parameter  of  injury.  For  the  tympanic 
membrane  and  the  upper  respiratory  tract,  the  critical  stress  calculated  using  the  mathematical 
model  differ  from  that  calculated  using  the  measured  data  by  an  amount  not  significantly  greater 
than  the  variation  of  stress  as  a  correlate  of  injury.  For  the  lung,  however,  the  variations  are  much 
larger  and  suggest  that  a  more  accurate  prediction  scheme  may  be  required. 
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1.  Parameter  Values  for  the  Three  Organ  Systems 


1.  INTRODUCTION 


The  repeated  passage  of  pressure  waves  past  a  given  point,  following  an  explosion  in  an 
enclosure,  results  in  the  long  and  complicated  pressure  time  histories  called  "complex  waves."  The 
waves  are  not  only  complex  because  of  their  temporal  variation,  but  their  effect  on  biological 
structures  depends  on  the  direction  of  each  wave  component.  The  ability  to  estimate  injury  under 
these  circumstances  represents  one  of  the  greatest  challenges  of  biomechanical  modeling. 

A  mathematical  model  for  the  reflections  of  blast  waves  within  an  enclosure  has  been 
constructed  based  on  the  concept  of  mirror-like  reflections  of  the  wave  front  at  each  solid  surface. 
The  model  accounts  for  the  position  and  orientation  of  each  surface  and  allows  for  the  presence  of 
openings,  such  as  doors  and  windows.  A  simplified  model  for  the  venting  of  explosion  product 
gases  is  also  included. 

The  pressure  loading  time  histories  are  combined  with  simplified  models  of  the  structural 
response  of  various  organ  systems  that  are  known  to  be  injured  by  air  blast.  The  combined  models 
not  only  allow  estimates  of  potential  injury  due  to  explosions  within  enclosures,  but  also  provide  a 
means  of  judging  the  accuracy  of  the  blast  calculations  in  terms  of  the  requirements  of  injury 
prediction. 

The  blast  field  predictions  are  compared  with  measurements  made  inside  an  Armored 
Personnel  Carrier  and  inside  a  field  bunker.  Many  qualitative  features  are  captured  but  the 
agreement  does  not  extent  to  individual  details.  A  comparison  of  the  peak  overpressures  and  the 
total  pressure  impulse  shows  that  in  most  cases  agreement  is  within  20-30%,  but  there  are  certain 
cases  for  which  the  disagreement  is  much  greater. 

The  measured  and  calculated  pressure  histories  are  used  in  the  simplified  biological 
response  models  and  a  comparison  made  of  the  critical  parameter  of  injury.  For  the  tympanic 
membrane  and  the  upper  respiratory  tract,  the  critical  stress  calculated  using  the  mathematical 
model  differ  from  that  calculated  using  the  measured  data  by  an  amount  not  significantly  greater 
than  the  variation  of  stress  as  a  correlate  of  injury.  For  the  lung,  however,  the  variations  are  much 
larger  and  suggest  that  a  more  accurate  prediction  scheme  may  be  required. 
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2.  COMPLEX  WAVE  MODEL 


The  first  step  in  COMPLX  is  to  use  the  method  of  images(3’'^)  to  generate  a  distribution  of 
changes  which  produces  the  same  wave  pattern  within  the  enclosure. 

This  algorithm  produces  reflected  Images  of  the  blast  source,  creating  N  primary  images 
(Fig.  1),  where  N  is  the  number  of  enclosure  walls.  The  primary  images  are  reflected  be  the  walls 
of  the  enclosure,  creating  N*(N  -  1)  secondary  images  (Fig.  2).  Images  do  not  reflect  off  the  wall 
that  created  them.  The  secondary  imager,  reflect  into  the  walls  of  the  enclosure  creating 
N*(N  —  1)2  tertiary  images  (Fig.  3).  This  is  an  infinite  process,  producing  reflecting  images  at 
increasing  distances  from  the  original  source.  The  model  carries  out  the  process  until  the  images 
reach  a  specific  distance  from  the  source.  Images  beyond  this  distance  do  not  contribute  to  the 
pressure-time  history  during  an  interval  equal  to  the  maximum  distance  divided  by  the  speed  of 
sound. 

First,  the  position  of  each  image  is  computed  as  a  function  of  the  position  of  its  source 
within  the  enclosure.  The  equations  used  for  computing  image  positions  are: 

Xi  =  Xj  -  (2  Djuff  Xjjorml) 

^i  “  *  (2  Djurf  Ynorml) 

^i  “  2^  ■  (2  Dsurf  ^norml) 

where: 

Xj  =  X  coordinate  of  Image 

Yj  =  Y  coordinate  of  Image 

Zj  =  Z  coordinate  of  Image 

Xj  =  X  coordinate  of  Image  source 

Yj  =  Y  coordinate  of  Image  source 

2^  =  Z  coordinate  of  Image  source 

^norml  ^  component  of  the  reflecting  wall’s  unit  vector 

^norml  =  ^  component  of  the  reflecting  wall’s  unit  vector 

Z^norml  ~  ^  component  of  the  reflecting  wall’s  unit  vector 

Djurf  =  distance  from  the  image  source  to  the  reflecting  wall. 
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Figure  1.  Primary  images  (Pj).  i  =  reflecting  wali. 


Figure  2.  Secondary  images  (Sjp.  i  =  primary  image  source; 
j  =  reflecting  wall. 
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Figure  3. 


Tertiary  images  (Tyi^.  ij  =  secondary  image  source; 
k  =  reflecting  wall. 
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Next,  the  pressure-time  history  is  computed  as  the  sum  of  all  blast  waves  that  can  be  viewed 
at  the  sensor  location  of  interest.  The  model  computes  the  path  from  the  sensor  back  to  the 
original  blast  via  each  image  and  its  source.  For  example,  a  tertiary  image’s  path  would  involve  the 
reflection  of  the  blast  off  three  walls  (Fig.  4).  COMPLX  checks  each  path  to  determine  if  it  can  be 
•Viewed"  from  the  sensor,  that  is,  that  the  path  is  not  obstructed  by  any  object  and  that  the  path 
intersects  the  reflecting  wall.  COMPLX  creates  a  list  of  all  paths  that  are  in  view  of  the  sensor. 

The  pressure-time  history  is  the  sum  of  the  individual  blast  waves,  each  of  which  is 
considered  to  be  a  Friedlander  wave.  The  characteristics  of  the  blast  waves  are  calculated  using  the 
normalized  blast  parameters  quoted  in  Chapter  6  of  Baker.(5) 


Pressure- time  history: 


k 

P(t)  -  ^  P^(t) 
i-1 


0  <  t  <  T 

max 


where: 

Tj^ax  given  by  the  maximum  specified  distance  an  image  can  be  from  the 
blast  divided  by  the  speed  of  sound. 

i  =  blast  wave  index 

k  =  number  of  blast  waves  seen  by  sensor 

[-b(t^/td)] 

Friedlander  wave:  P(t)-P  (l-t/t,)e 

s  a  d 

where: 

P(t)  =  pressure  at  time  t 

Pj  =  peak  pressure 

tg  =  time  measured  from  wave  arrival 

tj  =  positive  duration  of  wave 

b  =  exponential  parameter  of  Friedlander  wave 

The  parameters  Pj,  t^,  t<j,  and  b  can  be  determined  from  Baker  tables  when  the  normalized 
distance  from  the  sensor  is  known. 


Dj  -  d/D^ot 


3 


Figure  4.  Third  order  reflection.  The  distance,  ds,  between  the  sensor  (t)  and  the  tertiary 
image,  T4^2>  equal  to  the  path  of  the  blast  wave,  d4,  d3,  d2,  di . 

(ds  =  d4  +  d3  +  d2  +  di) 
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where: 


Dj  =  normalized  distance  from  source 
d  =  distance  of  blast  wave  path 
Dnot  =  (E/Po)'/3 
E  =  energy  released  by  charge 
Pq  =  ambient  pressure 

The  blast  parameters  t  and  t^j  are  determined  from  the  relations: 
t  =TaTo 
td  “  ^0 

where: 

Ta  =  normalized  time  of  arrival 
Tg  =  normalized  positive  duration 

Tq  =  Djiot/a 

a  =  ambient  speed  of  sound 

The  peak  pressure,  Pg,  is  determined  both  by  the  distance  to  and  strength  of  the  blast  and 
the  orientation  of  the  incident  wave  to  the  sensor. 

for  -  -  <  <  - 

Pg  =  [Pr  cos(e)  +  Pj  sin(0)]  Pq  R 
for  <?  >  I  f  I 

Pg  =  Pi  [1  +  0.25  cos(8)]  Po  R 

where: 


0  =  angle  between  the  individual  blast  waves  and  the  orientation  of  an  object  at  the  given 
location 


Pj  =  normalized  reflected  pressure 
Pj  =  normalized  incident  pressure 
Pq  =  ambient  pressure 
R  =  reflectivity  factor 

The  parameter  R  is  the  product  of  the  reflectivities  of  each  wall  the  wave  encounters  before 
reaching  the  sensor. 

COMPLX  adds  the  pressurization  of  the  enclosure  due  to  expanding  gases  of  the  blast  to 
the  pressure-time  history.  This  corresponds  to  the  "Quasi-Static  Overpressure  Response"^^)  '>nd 
"fill"!^).  The  pressurization  is  the  result  of  competition  between  heating  of  the  ambient  gases  in  the 
enclosure  by  the  explosive  and  the  loss  of  energy  from  the  flow  leaking  out  of  the  enclosure.  The 
deviation  for  computing  the  pressurization  of  the  enclosure  follows  as: 

Equation  of  state  for  the  entire  enclosure: 
pV  -  NRT 


E  -  E  +  NC  (T  -  T  ) 
O  V  o 


The  explosion  violently  produces  gas  which  adds  mass  and  energy  to  the  enclosure.  Both  factors 
increase  the  pressure  of  the  enclosure: 


dp  dE 

r '  e“ 


(1) 


The  explosion  adds  an  amount  of  energy,  E,  when  the  explosion  gases  completely  mix  with  the 
ambient  gases, 


V  (t) 

E(t)  -  E  +  AE  - - 

o  v„ 


If  the  mixing  proceeds  at  particle  speed,  Ug,  behind  the  blast,  then 


V  (t)  4/3  ’fCu  c)' 

E  e 


4/3 


u  t 
e 


-{rl 


where  the  volume,  V,  of  the  enclosure  is  represented  by  the  equivalent  sphere  and  Tg  is  the  time 
for  the  expanding  gases  to  fill  the  enclosure, 
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T  [t J 


for  t  <  T 


(2) 


o  e 


If  there  is  a  place  for  the  gases  to  leak  to  the  ambient,  then  the  flow  will  take  energy  out  of  the 
enclosure  at  a  rate 


M  _  zl 

IdtJ  V  ^out  out 

-'out 


The  flow  is  driven  by  the  energy  difference  between  the  outside  and  the  inside, 


u  -  kAp 
out 


so  that 


-  f— 1 

E  [dtj 


out 


k  p  A 
c  out 


^  a  i_  r^i 
p.  [p J 


Tq  depends  on  the  effectiveness  of  the  leak,  where: 


Tq  =  0  for  a  closed  enclosure 


Tq  =  oo  for  a  free  field. 


Combining  Eqs.  (1),  (2),  and  (3)  gives 


d  fApl  1  fdp^ 

dt  [p  J  T  [p  J 


o  o 


o  '  e 


Eq.  (4)  can  be  written  in  terms  of  dimensionless  variables. 


fApl  fAEl 

IpIJ  “  M 


/(t/t^) 


t/T  <  1 
e 


t/T  >  1 
e 


(3) 


(4) 


(5) 


where 


r  -  t/t 


to  yield 
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p'  + 


where 


o 


The  solution  is  then 


T  <  1  ' 

1  <  ^  . 


(«r)^  -  2(«r)  +  2  -  2  e  , 

> 

T  <  1 

V  -  2»  *  2  -  2  exp-'] 

► 

1  <  r 

Using  the  equations  of  state,  the  equation  for  the  expansion  pressure  becomes 

p  (t)  -  E  /V  f(0 
e  e  r 


(6) 


(7) 


where: 


Eg  =  energy  of  explosive 
Vj.  =  volume  of  enclosure 


The  final  step  in  computing  the  load  pressure-time  history  is  to  add  the  pressurization  of  the 
enclosure, 


p(t)  -  p(t)  +  p  (t) 
e 
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3.  ORGAN  RESPONSE  MODEL 


COMPLX  uses  the  load  pressure-time  history  determined  in  the  previous  section  as  the 
driving  force  for  the  Generalizable  ModelX^)  The  Generalizable  Model(2)  uses  the  damped 
harmonic  oscillator  equation  for  computing  the  structural  dynamics  of  the  body  due  to  a  blast 
wave.  Therefore  the  equation  of  motion  for  the  body  becomes: 


m 


dt^ 


r2m1  ^ 


+  kx  =  p(t) 


where: 

m  =  mass/area 

tg  =  characteristic  damping  time 

k  =  spring  constant/area 

X  =  displacement 

By  solving  for  the  displacement,  x,  and  the  velocity,  dx/dt,  COMPLX  can  compute  the 
delivered  stress  to  a  particular  organ.  The  model  uses  the  displacement  to  compute  the  stress 
delivered  to  the  tympanic  and  larynx  membranes  and  the  velocity  to  compute  the  stress  delivered 
to  the  lung. 

For  the  tympanic  and  larynx  membranes,  the  stress  is  computed  from: 

a  .  -  fkx 

tissue 

For  the  lung,  the  tissue  stress  is  assumed  to  be  proportional  to  the  parenchymal  pressure  at 
the  pleural  surface,  which  is  given  by: 

p  2m  dx 

p  t  dt  ’ 

c 

When  the  maximum  tissue  stress  within  an  organ  exceeds  a  critical  value,  injury  is  predicted 
to  occur.  Table  1  gives  the  parameter  values  for  the  organ  systems. 


Table  1.  Parameter  Values  for  the  Three  Organ  Systems 


Organ 

m 

(kgm/m2) 

^c 

(ms) 

k 

(kPa/mm) 

f 

®crit 

(MPa) 

Tympanic  membrane 

0.2 

10 

7.9 

133 

7.5-15 

Upper  respiratory  tract 

10 

100 

4 

5 

1 

Lung 

15 

10 

0 
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4.  COMPARISON  OF  COMPLEX  WAVE  MODEL 

WITH  DATA 


We  compare  the  predictions  of  COMPLX  with  two  sets  of  complex  blast  environment  data. 
The  first  set  was  measured  inside  an  APC  (Fig.  5).  The  second  was  measured  inside  a  field  bunker 
(Fig.  6).  The  digitized  field  data  is  plotted  as  a  pressure-time  history  curve  for  comparison  with  the 
same  curve  produced  by  the  model. 

Since  the  exact  conditions  of  the  test  are  not  known,  the  comparisons  are  only  approximate 
but  are  intended  to  give  the  reader  an  idea  of  COMPLX’s  ability  to  predict  the  characteristics  of  a 
complex  wave.  In  the  following  sequence  of  plots,  for  both  the  APC  and  the  bunker  data,  the 
parameters  of  the  model  are  set  up  such  that; 

1.  Sensor  positions  and  orientations  are  identical  to  those  of  the  field  blast. 

2.  The  orientation  function  is  the  same  as  described  in  the  second  section  of  this  report. 

3.  The  walls  have  100  percent  reflectivity. 

4.  The  enclosure  is  not  pressurized. 

The  APC  data  was  measured  from  a  series  of  blasts  in  which  the  explosive  charge  weights 
were  varied.  The  specific  charge  weights  were  57  gm,  113  gm,  227  gm,  and  454  gm.  The  pressure 
sensors  (Fig.  5)  used  to  measure  the  data  were  kept  in  the  same  location  for  the  entire  blast  series. 
The  N  and  S  pressure  sensors  faced  upwards  while  the  Lambdroid’s  four  sensors  faced  towards  the 
blast,  away  from  the  blast,  and  perpendicular  to  the  blast.  The  charges  were  all  detonated  in  the 
center  of  the  APC.  Comparison  plots  have  been  made  with  all  the  APC  data.  Figure  7  is  an 
example  from  this  set. 

The  bunker  data  was  measured  from  a  series  of  blasts  in  which  the  sensor’s  position  (Fig.  6) 
was  varied.  In  addition,  for  each  sensor  position  a  different  charge  weight  was  detonated,  227  gm 
and  454  gm.  The  sensor  configuration  (Fig.  6)  for  this  series  of  blasts  was  a  free  field  pointing 
upward,  Lambdroid  with  a  sensor  pointing  towards  the  blast,  a  sheep  with  a  skin  sensor  pointing 
towards  the  blast,  and  another  sheep  with  an  esophagus  sensor.  Comparison  plots  have  been  made 
with  all  Bunker  data.  Figure  8  is  an  example  from  this  set. 

We  have  taken  peak  pressure  and  total  impulse  from  each  of  the  comparison  data  sets  and 
plotted  field  data  versus  model  calculation  (Figs.  9  &  10).  This  creates  a  scatter  plot  such  that 
those  points  which  lie  on  the  45  degree  line  are  of  equivalent  values  for  field  data  and  model 
calculation.  The  points  that  don’t  lie  on  the  45  degree  line  are  off  by  a  percentage  difference. 
Scatter  plots  are  a  good  means  of  visualizing  the  correlation  of  a  large  data  set. 
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Figure  6.  Bunker  specification  and  sensor  location. 
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PRESSURE  HISTORY 


Rgure  7.  57  gm  C-4  blast  "S*  free  field  sensor. 
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5.  COMPARISON  OF  ORGAN  RESPONSE  USING 
CALCULATED  AND  MEASURED  PRESSURES 


Both  the  predicted  and  measured  pressure-time  history  curves  are  used  as  input  to  the 
Generalizable  ModeK^)  to  produce  a  sequence  of  three  plots.  The  first  plot  graphs  the  maximum 
tympanic  membrane  stress  as  a  function  of  time.  The  second  plot  graphs  the  maximum  larynx 
stress  as  a  function  of  time.  The  third  plot  graphs  the  lung’s  overpressure  as  a  function  of  time. 
These  quantities  are  used  in  the  prediction  of  blast  injury. 

Figures  1 1  and  12  are  the  predicted  organ  response  where  the  driving  forces  are  plotted  in 
Figures  8  and  9,  respectively.  Scatter  plots  have  been  made  for  the  maximum  tympanic  stress  (Fig. 
13),  the  maximum  URT  stress  (Fig.  14),  and  the  work  done  on  the  lung  (Fig.  15).  A  case-by-case 
description  may  be  found  in  a  supplementary  report. 


Field  Data 


Model 


Figure  11.  57  gm  C-4  blast  "S'  free  field  sensor. 


Model 


Figure  12.  2/1 3/87  blast  Lambdroid  face-on  sensor. 
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Figure  13.  Tympanic  stress. 
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